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Abstract: The robustness of the large mixing angle (LMA) oscillation (OSC) interpre- 
tation of the solar neutrino data is considered in a more general framework where non- 
Qh, standard neutrino interactions (NSI) are present. Such interactions may be regarded as 

Q_i| a generic feature of models of neutrino mass. The 766.3 ton-yr data sample of the Kam- 

LAND collaboration are included in the analysis, paying attention to the background from 
the reaction 13 C(a, n) 16 0. Similarly, the latest solar neutrino fluxes from the SNO collabo- 



ration are included. In addition to the solution which holds in the absence of NSI (LMA-I) 
there is a "dark-side" solution (LMA-D) with sin 2 9 S0L = 0.70, essentially degenerate with 
the former, and another light-side solution (LMA-0) allowed only at 97% CL. More precise 
KamLAND reactor measurements will not resolve the ambiguity in the determination of 
the solar neutrino mixing angle 6 S ol, as they are expected to constrain mainly Amg 0L . We 
comment on the complementary role of atmospheric, laboratory (e. g. CHARM) and fu- 
ture solar neutrino experiments in lifting the degeneracy between the LMA-I and LMA-D 
solutions. In particular, we show how the LMA-D solution induced by the simplest NSI 
between neutrinos and down-type-quarks-only is in conflict with the combination of current 
atmospheric data and data of the CHARM experiment. We also mention that establishing 
the issue of robustness of the oscillation picture in the most general case will require further 
experiments, such as those involving low energy solar neutrinos. 
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1. Introduction 

The very first data of the KamLAND collaboration [1] have been enough to isolate neutrino 
oscillations as the correct mechanism explaining the solar neutrino problem [2], indicating 
also that large mixing angle (LMA) was the right solution. The 766.3 ton-yr KamLAND 
data sample strengthens the validity of the LMA oscillation interpretation of the data [3]. 

With neutrino experiments now entering the precision age [4], the determination of 
neutrino parameters and their theoretical impact have become one of the main goals in 
astroparticle and high energy physics [5]. Now the main efforts should be devoted to the pre- 
cision determination of the oscillation parameters and to test for sub-leading non-oscillation 
effects such as spin-flavour conversions [6, 7] or non-standard neutrino interactions (NSI, 
for short) [8]. 

A quantitative analysis of neutrino oscillations reveals that the interpretation is rela- 
tively robust even taking into account the possibility of solar density fluctuations in the 
solar radiative zone [9], that might arise from magnetic fields effects [10], currently uncon- 
strained by helioseismology. The robustness of neutrino oscillations in the presence of spin- 
flavour conversions induced by non- vanishing neutrino transition magnetic moments [11] 
follows from the stringent limit on anti-neutrinos from the Sun by the KamLAND collab- 
oration [12] . 



1 This does not hold for the Dirac case, but here the theoretical expectations for magnetic moments are 
typically much lower. 



Here we focus on the case of neutrinos endowed with non-standard interactions. These 
are a natural outcome of many neutrino mass models [13] and can be of two types: flavour- 
changing (FC) and non- universal (NU). 

Seesaw-type models leads to a non-trivial structure of the lepton mixing matrix char- 
acterizing the charged and neutral current weak interactions [14]. This leads to gauge- 
induced NSI which may violate lepton flavor and CP even with massless neutrinos [15-19]. 
Alternatively, non-standard neutrino interactions may also arise in models where neutrino 
masses are "calculable" from radiative corrections [20, 21]. Finally, in some supersymmet- 
ric unified models, the strength of non-standard neutrino interactions may be a calculable 
renormalization effect [22]. 

How sizable are non-standard interactions will be a model-dependent issue. In some 
models NSI strengths are too small to be relevant for neutrino propagation, because they 
are suppressed by some large scale and/or restricted by limits on neutrino masses. However, 
this need not be the case, and there are interesting models where moderate strength NSI 
remain in the limit of light (or even massless) neutrinos [15-19]. Such may occur even in 
the context of fully unified models like SO(10) [23]. 

Non-standard interactions may in principle affect neutrino propagation properties in 
matter as well as detection cross sections [2]. Thus their existence can modify the solar 
neutrino signal observed at experiments. They may be parametrized with the effective 
low-energy four-fermion operator: 

Cnsi = -e f a pV2G F (va^Lvp) (f-fPf) , (1.1) 

where P = L, R and / is a first generation fermion: e, u, d. The coefficients sr a denote the 
strength of the NSI between the neutrinos of flavours a and (5 and the P-handed component 
of the fermion /. In the present work, for definiteness, we take for / the down- type quark. 
However, one can also consider the presence of NSI with electrons and up and down quarks 
simultaneously. Current limits and perspectives in the case of NSI with electrons have been 
reported in the literature [24]. 

While strong constraints exist from v^ interactions with a down-type quark (e™T < 
10_3 , £ % ^ 10_3 - 10_2 ) from CHARM and NuTeV [25], the constraints for all other 
NSI couplings, including those involved in solar neutrino physics, are rather loose [25, 26]. 
Therefore, in our analysis we consider e~T = and we concentrate our efforts in the rest 
of NSI parameters. 

For our solar neutrino analysis, we will consider the simplest approximate two-neutrino 
picture, which is justified in view of the stringent limits on #13 [5] that follow mainly from 
reactor neutrino experiments [27]. 

The Hamiltonian describing solar neutrino evolution in the presence of NSI contains, 
in addition to the standard oscillations term 

(-^co S 26 + V2G F N e 4^sin2#\ 

lv ^ sin 20 ^ cos 20 J ^ 



a term H^si, accounting for an effective potential induced by the NSI with matter, which 
may be written as: 

H NS i = V2G F N d (° £ £ £l ) ■ (1.3) 

Here e and e' are two effective parameters that, according to the current bounds 

f P 

discussed above (e afJl ~ 0) , are related with the vectorial couplings which affect the neutrino 
propagation by: 

e = - sin 9 23 ef T e' = sin 2 9 23 ef T - ef e (f .4) 

The quantity N d in Eq. ( |f.3| ) is the number density of the down-type quark along the 
neutrino path. In the more general case, the effective couplings e and e' will contain 
contributions from the three fundamental fermions and NSI effects would be important 
not only in neutrino propagation but also in the detection process. 

It is important to note that the neutrino evolution inside the Sun and the Earth is 
sensitive only to the vector component of the NSI, e^ = s& + e^f. The effect of the axial 
coupling will be discussed in detail in section 4. 

Before introducing our numerical analysis of the solar neutrino data in the next sec- 
tion, it is worth discussing the analytical formulas for neutrino survival probability in the 
constant matter density case, in order to have a better understanding of the results that 
will be shown in the next section. Recall first that in two-neutrino oscillations [28, 29] 
one can, without loss of generality, restrict the variation of the mixing angle 9 only to the 
range [0, ^] and still cover the whole physical space 2 . In the adiabatic regime the survival 
probability can be approximated by Parke's formula [30] 

P(zy e -^ e ) = -[l + cos20cos20 m ], (1.5) 

where 9 m is the effective mixing angle at the neutrino production point inside the sun. In 
the absence of non-standard neutrino-matter interactions the mixing angle in matter may 
be obtained from the expression 

Am 2 cos29-2V2 EG F N e 
cos2# m = — = , (1.6) 



[Am 2 cos 2(9 -2^/2 EG F N e ) + (Am 2 sm29y 

In order to explain the deficit of solar neutrinos observed at the detectors, the neutrino 
survival probability should satisfy: P < 0.5. According to Eqs. Ql.5| ) and (|1.6|), this 
requirement is not satisfied for cos 29 < 0, so that only vacuum mixing angles in the "light 
side" (0 < 9 < ?) can give solution to the solar neutrino problem. Indeed this is confirmed 
by the results shown in Fig. ffl. 

Lets now turn to the case where non-standard interactions are present, in addition to 
oscillations. Within such generalized picture (OSC+NSI), Eq. (|1.6|) is modified to 



Am 2 cos29-2V2 EG F (N e -e'N d ) 
cos20 m = —^ , (1.7) 



2 Alternatively one can restrict the angle to the range < 9 < j if one includes a separate region with 
Am 2 < 0. As discussed in [28] it is more natural to use < 6 < ^ with a fixed sign of Am 2 . 



where 



[Am 2 ] 2 matter = [Am 2 cos 29 - 2^2 EG F (N e - e'N d )\ * 



+ 



Am 2 sin 29 + iV2eEG F N d 



(1-8) 



Thanks to the presence of the non-universal coupling e' one can obtain P < 0.5 even 

for cos 29 < as long as e' > F e ™ |cos — '-, This makes it possible to explain the 

solar neutrino data for values of the vacuum mixing angle in the dark side, for large enough 
values of e' . As we will see below this possibility will lead to the appearance of the LMA-D 
solution with 9 > j and thus to the ambiguous determination of the solar mixing angle. 

2. The fit 

Here we reanalyse the robustness of the oscillation interpretation of the solar neutrino data 
in the presence of non-standard interactions. We include the recent SNO data [31] as well 
as the 766.3 ton-yr data sample from the KamLAND collaboration [3], taking into account 
the background from the reaction 13 C(a,n) ie O. In order to do this we first calibrate 
with the results obtained in the pure oscillation case for the KamLAND-only, solar-only 
and combined data samples. For the KamLAND analysis we use a Poisson statistics as 
described in [32]. Our best fit point is located at sin 2 26> SO l = 0.84 and Am 2 0L = 7.9 x 10 -5 
eV 2 , in good agreement with the results of Ref. [3] . For the solar data we include the rates 
for the Chlorine, Gallex/GNO, SAGE, as well as the Super-Kamiokande spectrum, SNO 
day/night spectrum and SNO salt data, with a total of 84 observables. We adopt the pull 
method [33] to fit the data using the most recent BS05 Standard Solar Model [34]. We 
perform a complete analysis of the solar neutrino data using a numerical computation for 
the survival probabilities both in the light as well as in the dark side of the mixing angle, 
for values of Am 2 0L in the range of 10 -6 to 10 -3 eV 2 and running also both e and e' at the 
same time in the range [—1, 1]. Our results for e = e' = are shown in Fig. ffl. The best 
fit point for this global analysis is given by sin 9 S0L = 0.29 and Am 2 0L = 8.1 x 10 -5 eV 2 . 
This is in excellent agreement with the results obtained in [5] for the solar case. Reassured 
by this calibration we now turn to the generalized OSC+NSI picture. Our results for this 
case are shown in Fig. || and Table [l]. One sees that, in the light side, we obtain a region of 
allowed oscillation parameters larger than in the pure oscillation case, but more restricted 
than those obtained in previous OSC+NSI analysis of Refs. [35, 36] due to the effect of 
the recent KamLAND data, visible mainly in Am 2 OL . The table gives the parameter best 
fit values for the OSC and OSC+NSI fits. For the OSC+NSI analysis the best fit occurs 
for e = 0.0 and e' = —0.05. Clearly the quality of the fit obtained with and without 
NSI is comparable, as seen from the x 2 values given in the last column of the table. The 
most remarkable result is, however, the appearance of an additional solution in the dark 
side region, which can be qualitatively understood from the discussion given at the end of 
Sect. [l]. This LMA-D solution has sin 2 9 S0L = 0.70 and the same Am, 2 0L value as the LMA-I 
solution and is significantly better than the LMA-0 OSC+NSI solution of Ref. [35, 36], as 
shown in the table. On the other hand, it is nearly degenerate with the LMA-I solution, 
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Figure 1: 90%, 95%, 99% and 99.73% C.L. allowed regions of the neutrino oscillation parameters 
from the analysis of the latest solar data (hollow lines, left panel), and latest KamLAND data 
(colored regions, left panel) and from the combined analysis (right panel). 
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Figure 2: Allowed regions for the generalized OSC + NSI case, determined from the latest data: 
left panel corresponds to a solar only analysis, while the right panel corresponds to the combined 
solar+KamLAND analysis. 



as seen by the \ 2 value. This solution is characterized by e' = 0.90, although lower values 
~ 0.75 are allowed at 3a. Although embarrassingly large, one sees that such large NSI 
strength values are perfectly compatible with all existing solar and reactor neutrino data, 
including the small values of the neutrino masses indicated by current oscillation data. This 
opens a potentially physics challenge for upcoming low energy solar neutrino experiments, 
such as Borexino. Note that large NSI values could affect also solar neutrino detection, as 
considered in [37]. In what follows we give a discussion of the role of other experiments in 
probing neutrino properties at the level implied by the above LMA-D solution. 
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Table 1: Best fit solar neutrino oscillation points with and without non-standard neutrino inter- 
actions. 




Figure 3: Constraining NSI parameters: dependence of A\ 2 with respect to e and e' , illustrating 
the current limits. 

3. Constraints on NSI: present and future 

As we just saw there are constraints on non-standard neutrino interaction strength param- 
eters that follow from current solar and KamLAND data. The existence of NSI could also 
potentially affect neutrino-nucleon scattering and there are laboratory data that poten- 
tially constrain their allowed strength. Moreover, one must check restrictions that follow 
from atmospheric data. Here we discuss their complementarity. 



3.1 Solar and KamLAND 

We can derive limits on NSI parameters from solar and KamLAND data by displaying 
our x 2 as a function of the NSI parameters e or e' and marginalizing with respect to 
the remaining three parameters. Figure y gives the A\ 2 profiles with respect to e and 
s'. From here one can determine the corresponding constraints on e and e' . We can see 
that at 90% C.L. -0.93 < e < 0.30 while for e' the only forbidden region is [0.20,0.78]. 
Note that our limits on e' are weaker than those of Ref. [36], which apply only to the 
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Figure 4: Predicted neutrino survival probability for low-energy neutrinos (left) and boron 
neutrinos (right) at the best fit points of LMA-I, LMA-D and LMA-0. 

restricted case where e = 0. We see that the limits on the strength of non-standard neutrino 
interactions are still very poor. The dashed lines in Fig. || denote the ultimate reach of this 
method of constraining NSI parameters (through their effect in solar neutrino propagation), 
namely they correspond to the case where solar neutrino oscillation parameters Amg 0L and 
6 S0L are determined with infinite precision. One sees that in this ideal case the allowed 
range narrows down mainly for negative NSI parameter values. We conclude that there 
is substantial room still left for sub-leading non-standard neutrinos conversions in matter 
and, moreover, that the determination of solar neutrino oscillation parameters, especially 
the solar mixing angle, is currently ambiguous. It is unlikely that more precise reactor 
measurements by KamLAND will resolve this mixing angle ambiguity, as they are expected 
to constrain mainly Am,g 0L . 

In Fig. H] we present the predicted neutrino survival probabilities versus energy, from 
the region of pp neutrinos up to the high energy solar neutrinos, for the three best-fit 
points of the allowed regions found above. One sees that the solutions predict different 
rates for the low energy neutrinos (e.g. pp and pep), so that future low energy solar 
neutrino experiments may have a hope of disentangling these solutions. Similarly, in the 
region of boron neutrinos our LMA-D solution also predicts a distortion in the spectrum 
that might be detectable at future water Cerenkov experiments such as UNO or Hyper- 
K [38], given the high statistics expected. With good luck such high statistics experiments 
may have a window of opportunity. 



3.2 Laboratory experiments 

The laboratory bounds on the neutrino non-standard interactions with down-type quarks 
can be summarized as |e^f | < 0.5, \eff\ < 6, |e^| < 1.1, -0.6 < e*f < 0.5, -0.3 < e*£ < 
0.3 3 , see e. g. Ref [25]. Here we are interested in vector-like NSI couplings. For the case 



3 There is a second branch 0.6 < eff" < 1.1 which should be added to the ranges given in Ref [25]. 




Figure 5: Consistency between the e' coupling required for our LMA-D solution (shaded band) 
and the regions allowed by atmospheric data in the analytic approximation of Ref. [40] for e^. = 0.21 
(solid lines). The laboratory constraints are also shown (dashed lines). See the text for a detailed 
explanation. 
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of efg 7 , these bounds can be translated to 
wider range. However, we stress that these bounds have been obtained assuming that only 
one parameter is effective at a time. Relaxing this assumption opens more freedom. This 
is why we have chosen to indicate them by the dashed lines in Fig. IEL Assuming maximal 
mixing in the 2-3 sector in Eq. (|l.4| ), one has 
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^From this one can see explicitly that, even taking the above constraints at face value, they 
still leave room for our degenerate dark-side solution with e' = 0.90. 



3.3 Atmospheric data 

Concerning the atmospheric neutrino data, it is known that a large NSI strengths can 
originate a suppression of the neutrino oscillation amplitude. This has indeed been used in 
a two-neutrino analysis [39] in order to obtain relatively strong bounds on the NSI strength. 
However, in a 3-neutrino analysis of atmospheric data [40] it has been explicitly shown that 
large NSI strengths are not excluded. In particular, these authors have found two specific 
scenarios where somewhat large NSI strengths can fit well the experimental data, because 
their effect will be indistinguishable from the standard oscillation case, at least at high 
and low energies. Adapting their definitions to our notation, and using their analytical 
description 4 , we obtain the two branches indicated in Fig. |g. One sees that the shaded 
band corresponding to our LMA-D solution at 90% C.L. (with e^. = 0.21) intersects these 
branches in two disjoint regions, suggesting that, indeed, the NSI couplings required by 



4 This analytical comparison holds only at high energies, a more complete check would require a numerical 
analysis to see the effect of intermediate energies. 



the LMA-D solution are compatible with the atmospheric neutrino data. However, in a 
more complete numerical analysis of atmospheric neutrino data [41], it has been shown 
that values of e^. in the right region are not allowed by atmospheric data: only the left 
disjoint region is compatible with atmospheric neutrino data. As indicated by the dashed 
lines in Fig. [5] one can see that e^T values in this region lie outside the range allowed by 
current laboratory data. This leads us to conclude that the LMA-D solution induced by 
the simplest non-standard interactions of neutrinos with only down-type quarks is ruled 
out by its incompatibility with atmospheric and laboratory data. However, one can verify 
that for the general case where neutrinos have other NSI couplings one can reconcile the 
above laboratory bounds with the parameters required by the LMA-D solution. 

Finally, we comment on the magnitude of flavor-changing NSI of neutrinos. First 
note that direct constraints on the magnitude of these couplings do not exist. The only 
"bounds" mentioned in the literature are obtained from charged lepton flavor violating 
processes. While the constraints are very restrictive, they are theoretically fragile, to the 
extent that they rely on the assumption of weak SU(2) symmetry, and can therefore be 
avoided if one allows for SU(2) breaking. Note, however, that the magnitude of the FC 
neutrino NSI required for our dark-side solution is quite small. Futuristic proposals for 
improving these constraints with coherent neutrino scattering off nuclei have been already 
discussed [42]. 

4. A comment on axial NSI couplings 

Before we conclude, let us mention that, up to now we have only considered the effects of 
the NSI on the neutrino propagation through the Earth and solar interior. These effects 
appear as a result of the vectorial couplings of neutrinos with down-type quarks. An axial 
component of the NSI coupling could give rise to a non-standard contribution to the NC 
cross section detection at the SNO experiment. As already noted in Ref.[25], the SNO-NC 
signal will be modified as: 

<t>NC~f B (l + 2e A ) (4.1) 

where 

SA = - Y, ( P ^N C 4i ( 4 - 2 ) 

with e d J^ = £^ — e^? being the couplings which enter into the effective Lagrangian. Thus, 



ea is independent of the effective couplings e and e' defined in Eq. (1.4). In the analysis 
performed so far we have assumed ea = 0. This assumption is well justified due to the 
good agreement between the SNO NC measurement and the SSM prediction for the boron 
flux. However, we now relax this assumption and include the effect of the new parameter 
ea in our analysis. 

The results obtained in a generalized analysis which takes into account the presence 
of an non-zero axial component of the NSI (5 parameters instead of 4) are summarized 
in Fig. ||. One sees that the neutrino data clearly prefers ea ~ 0, in agreement with our 
previous approximation, in Sec. 0. 
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Figure 6: Results of the analysis including the axial component of the NSI. In the left panel we 
compare the allowed regions at 90% and 3<r obtained with (color regions) and without (lines) ea- 
In the right panel we show the dependence of A% 2 with respect to the axial coupling ea 



5. Conclusions 

In short, we have reanalysed the status of the LMA oscillation interpretation of the so- 
lar neutrino data in a more general framework where non-standard neutrino interactions 
are present. We have seen that combining the solar neutrino data, including the latest 
SNO fluxes of the salt phase with the full KamLAND data sample still leaves room for 
a degenerate determination of solar neutrino oscillation parameters. To this extent the 
solar neutrino oscillation parameters extracted from the experiments may be regarded as 
non-robust. In addition to the lower LMA-0 solution, we have found a LMA-D solution 
characterized by values of the solar mixing angle larger than 7r/4. This solution requires 
large non-universal neutrino interactions on down- type quarks. While the LMA-0 solution 
is already disfavored, and will soon be in conflict with further data, e.g. future KamLAND 
reactor data, the degeneracy implied by LMA-D solution will not be resolved by more 
precise KamLAND reactor measurements. This shows that the determination of solar neu- 
trino parameters only from solar and KamLAND data is not fully robust. It is crucial to 
consider other data samples, such as atmospheric and laboratory data, since these bring 
complementary information. In the present case they allow one to rule out the LMA-D 
solution induced by the simplest NSI between neutrinos and down-type-quarks-only, given 
the large values of the non-universal NSI couplings required by that solution. It is therefore 
important to perform similar analyses for the more general case of non-standard interac- 
tions involving electrons and/or up-type quarks. Only in such scenario (NSI with u-type, 
d-type and electrons) we can confidently establish the robustness of the oscillation inter- 
pretation. Further experiments, like low-energy solar neutrino experiments are therefore 
required in order to clear up the situation. 
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